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Abstract 

A search is performed for production of short-lived particles in e'''e~ XY, with X — > 77 
and Y — > ff, for scalar X and scalar or vector Y. Model-independent limits in the range of 
25-60 femtobarns are presented on cr{e~^e~ XY) x B(X 77) x B{Y ff) for centre-of-mass 
energies in the range 205 — 207 GeV. The data from all LEP centre-of-mass energies 88 — 209 GeV 
are also interpreted in the context of fermiophobic Higgs boson models, for which a lower mass 
limit of 105.5 GeV is obtained for a "benchmark" fermiophobic Higgs boson. 
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1 Introduction 



This paper presents the results of two types of search for the production of a di-photon system 
recoihng from another massive scalar or vector object. The searches are sensitive to the processes 
e"'"e~ — > XY, with X — > 77 and Y — > ff, where ff is a hadronic system (jets), a pair of charged 
leptons, or neutrinos resulting in missing energy. In the general search mode X must be a scalar, 
Y can be any scalar or vector particle of any mass, and both particles must be short-lived so 
that they decay close to the interaction point. The other search mode is referred to as the h^Z^ 
search; it requires Y to be a boson and is applied to data taken at all energies. The data 
used for these searches were recorded by the OPAL detector at centre-of-mass energies (-E'cm) 
88 — 209 GeV, the entire energy range achieved at LEP. 

These searches are largely motivated by "fermiophobic" scenarios where one of the Higgs 
bosons decays primarily into a boson pair. In the fermiophobic interpretation, Y would be a 

and X a Higgs boson decaying into two photons. Indeed, the Higgs boson predicted in the 
Standard Model decays into two photons via a quark- or W-boson loop 0, but with a rate too 
low for observation of the process at LEP luminosities. Processes e'''e~ — > h'^Z^ — > 77ff at near- 
Standard Model production rate and having large di-photon branching ratios have been predicted 
in a number of alternative theories here h*^ refers to the lightest neutral boson where 

extended Higgs sector models are discussed. A particularly natural situation for fermiophobic 
Higgs bosons occurs in two Higgs doublet models (2HDM) of "Type-I", where one Higgs 
doublet couples only to bosons. Because there are different fermiophobic models, it is not 
possible to present search results for the entire parameter space of the various theories. In the 
present paper a benchmark fermiophobic model is defined as having Standard Model production 
strength and a Higgs boson di-photon branching fraction calculated by turning off the fermion 
couplings to the Higgs boson in the Standard Model. 

The OPAL Collaboration has presented searches similar to those reported here for LEP 
energies up to Ecn, = 189 GeV g-llll; this paper extends those searches with the addition of 
data taken at -Ecm = 192—209 GeV. Fermiophobic Higgs boson searches have also been presented 
by other LEP collaborations [p!^14| and by hadron collider experiments [15, If]. To date, no 
evidence of a fermiophobic Higgs boson has been seen. 



2 Data, Simulated Backgrounds and Signals 



The data used in this analysis were recorded using the OPAL detector [17| at LEP. The 1999 
data consisted of 217.0 ± 0.7 pb"^ collected at Ecm = 192 - 202 GeV. The 2000 LEP data 
consisted of 211.1 it 0.8 pb^^ collected at Ecm = 200 — 209 GeV, with the majority of the data 
taken at 205 and 207 GeV. The data sets are summarized in Table |l|. 

The backgrounds from Standard Model processes were modelled using Monte Carlo simula- 
tions at ^ = 192, 196, 202, and 206 GeV for the 1999 and 2000 data. Simulated events were 
processed using the full OPAL detector Monte Carlo ||l^ and analysed in the same manner as 
the data. The full-detector simulations were reweighted for the ^/s distribution of the data using 
the Monte Carlo generators. 

The dominant background to this search arises from the emission of two energetic initial state 
radiation (ISR) photons. This process was simulated using the KK2f/CEEX generator with 



hadronisation and fragmentation by PYTHIA 6.125 |2C]. The CEEX modelling of ISR employs 
full second-order QED corrections to the matrix element, and applies coherent exponentiation of 
the QED corrections from interference between ISR and final-state radiation. Four-fermion pro- 



cesses were modelled using the grc4f |21| and KORALW [22| generators. Two-lepton final states 
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were simulated using BHWIDE |23|], TEEGG || and KORALZ |25|]. The NUNUGPV 
program was used to generate events of the type e+e" — > 1/2^77(7). The process e+e" — > 77 was 
simulated using the RADCOR generator [27|. Tau lepton decays were modelled using Tauola 
2.4 



The process e~''e - ^ h^zo, hO — > 77 was simulated for each Z° decay channel using the 



HZHA3 generator |29|. For the general search, which is applied to the data taken in 2000 only, 
the role of the was replaced by scalar or vector particles having masses from 10 — 200 GeV. 
Efficiencies for signals were estimated by generating Monte Carlo for both scalar and vector 
signals in mass steps of 5 GeV; the scalar and vector efficiencies agreed within systematic errors 
and therefore are not treated as separate cases. 



3 Event Selection 



The analysis described in the following is identical to the one used in the paper for OPAL data 
taken at 189 GeV p. Slightly different analysis cuts were used on the lower energy data sets, 
as described in the earlier publications [P-|fl|- 

Events were selected if there were at least two photons recoiling from some other system 
decaying into one of the following three topologies: 



(A) a qq pair ("Hadronic Channel"), or 

(B) one or two charged leptons ("Leptonic Channel"), or 

(C) a z^z? pair ("Missing Energy Channel"). 



Photons were identified as clusters in the electromagnetic calorimeter (EC) which were not 
associated with tracks if the lateral spread of the clusters satisfied the criteria described in 



reference |10|. The efficiencies were increased by approximately 10 — 20% by including photon 



conversions into e"'"e pairs using the methods described in reference 

The dominant background to the searches arises from ISR producing mostly low-energy 
photons along the beam direction. Therefore, we required the two highest-energy photons in 
the event to satisfy the following: 



(Gl) The two photon candidates were required to be in the fiducial region | cos(0^)| < 0.875, 
where the polar angle is the angle of the photon with respect to the e~ beam direction. 

(G2) The highest-energy photon was required to have -E'71 /-Ebeam > 0.10 and the second- highest- 
energy photon was required to have -E72 /-E-beam > 0.05. 

(G3) The sum of track momenta and extra electromagnetic cluster energies in a 15 degree cone 
about the photons had to be less than 2 GeV. 



The remaining cuts depend on the particular recoil topologies. In order to assess the 
background modelling, the photon cuts are not applied until after preselection cuts for the 
three final state topologies. For all topologies, tracks and EC clusters that are not associated to 
tracks are required to satisfy the criteria defined in reference . The criteria for the definition 
of tracks and EC clusters in the OPAL detector are described in reference [31|. 
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3.1 Hadronic Channel 



The hadronic channel is characterised by two photons recoihng against a hadronic system. 
Candidate events were required to satisfy the foUowing criteria: 

(Al) The standard OPAL hadronic event preselection in Ref. i?vis > 0.5; |S p2^\ < 

0.6£'beam; and at least two electromagnetic clusters with E^/E'beam > 0.05. The quantities 

E'vis and pVis are the scalar and vector sums of track momenta, unassociated EC and 

E 

unassociated hadron calorimeter cluster energies, and i?vis = '^^^ visible momentum 

along the beam direction, obtained from the sum of all tracks and unassociated clusters, 
is denoted by |S Pz'^^l- 
(A2) The photon pair criteria Gl— G3. 

(A3) Photon isolation: both photon candidates were required to satisfy px, jet-7 > 5 GeV, 
where pT, jet-7 photon momentum transverse to the axis of the closest jet out of 

two jets formed with the Durham scheme (excluding the photon pair). 

(A4) Photon energy balance: {Eyi - E^2)/Eo < 0.5, where Eo = {s - M^)/{2^/s) would be the 
energy of a single photon recoiling from the TP . This cut discriminates against ISR photon 
pairs. 

(A5) The recoil mass from the di-photon system, Mrecoib is required to be consistent with the 
TP: |Mrecoii — -^z| < 20 GcV (not used in the general search mode). 



3.2 Charged Lepton Channel 

This channel searches for events in the 77^^^" final state. Events having only one well-identified 
lepton are accepted to avoid efficiency loss for lepton tracks at low polar angles. The lepton 
tracks are treated as jets to include tau lepton final state topologies. Leptonic channel candidates 
were required to satisfy the following criteria: 



(Bl) The standard OPAL low multiphcity preselection of Ref. IJ; i?vis > 0.2; |S p™| < 
0.8£'beami number of EC clusters not associated with tracks < 10; number of tracks A^t 
satisfies 1 < At < 7; at least two electromagnetic clusters with £'/Sbeam > 0.05. 

(B2) The photon pair criteria Gl— G3. 

(B3) For events having only one track and a converted photon, the EC cluster associated with 

the track must not also be associated with the conversion. 
(B4) For events having two or more tracks, the event is forced to have two jets within the 

Durham scheme, excluding the identified di-photon candidate, and both jets are required 

to have energies above 3 GeV. 
(B5) iMrecoii " Afz| < 20 GcV (not used in the general search mode). 



3.3 Missing Energy Channel 



The missing energy channel is characterised by two photons and no other significant detector 
activity. Candidates in the missing energy channel were required to satisfy the following criteria: 



(Gl) The standard OPAL low multiplicity preselection of Ref. [34|; the vetoes in Ref. 

against cosmic ray and beam- wall/beam-gas backgrounds; number of EC clusters not 
associated with tracks < 4; number of tracks < 3; |S p^^\ < 0.8£'beamj and at least two 
electromagnetic clusters with i?/£'beam > 0.05. 
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(C2) The photon pair criteria Gl— G3. 

(C3) pxini) > O.OSii^beam where Pt{ii) is the transverse momentum of the di-photon system; 
the angle between the two photons in the plane transverse to the beam axis: |0-y-y — 180° | > 
2.5°; the polar angle of the momentum of the di-photon system: | cos0^^| < 0.966. 

(C4) No track candidates other than those associated with an identified photon conversion. 

(C5) Veto on unassociated calorimeter energy: the energy observed in the EC not associated 
with the two photons is required to be less than 3 GeV. 

(C6) iMfecoii — -^zl < 20 GeV (not used in the general search mode). 



4 Results 

For the 1999 and 2000 data, the numbers of events passing the cuts are listed for the three recoil 
topologies in Table |. There are no events in which more than one photon pair satisfying the 
cuts was found. The numbers of candidates passing cuts are generally in good agreement with 
the expected numbers of Standard Model backgrounds; this was also the case in earlier OPAL 



publications for the lower Ecm |pHll|- The one noteworthy discrepancy is for cut CI in the 
missing energy channel. In this channel there is a large background from Bhabha electrons lost 
in the beampipe. The ISR photons for this background have a steeply rising population in the 
forward direction. Cut CI is made before the cut on polar angle, and therefore a discrepancy 
arises because of the steep angular distribution and the inadequate modelling of material in the 
very low polar angle regions. 

Combining both the 1999 and 2000 data in the three topologies, 112 candidates pass the 
general cuts compared to 118.3ib7.9 expected background, and 42 candidates pass the h°ZOcuts 
compared to 51.9ib2.9 expected background. 



4.1 Systematic Errors 

The uncertainty on the modelling of ISR is the most important component of the systematic error 
because of the irreducible background arising from this process. This uncertainty is estimated 
from the comparison of data with the Standard Model background simulation for events passing 
cuts A2, B2, or C2. The shapes of the distributions for E^i and E^2 are modelled well by the 
simulations. The simulations also reproduce well the number of events observed in the three 
channels combined. For the 1999 and 2000 runs the simulations predict 3.7% and 5.0% fewer 
events than observed, respectively. The statistical error on the 1999 data is 4%, and similarly 
for the 2000 data. Modelling of photon conversions has an uncertainty of approximately 1%. 
Uncertainties on the integrated luminosities of the data sets are negligible compared to the other 
uncertainties. Combined, these error sources result in a total background uncertainty estimate 
of 10%. The experimental results which follow are not very sensitive to this number. 

The dominant systematic uncertainty for the signal acceptances arises from the photon 
detection efficiency, primarily due to the simulation of the photon isolation criterion G3 [ pTj] , 
and is estimated to be 3%. The uncertainty from Monte Carlo statistics is typically better than 
4%. A systematic error on the photon energy scale is estimated by comparing the fitted single- 
photon ISR energy peak to the expected value based on the precisely known beam energy and 
TP mass. For this study a sample of single-photon events was generated and compared to the 
data for photon energies above 5 GeV and polar angles greater than 25 degrees. This leads to a 
systematic uncertainty on the di-photon mass of 0.35 GeV at a mass of 100 GeV. The resolution 
on the di-photon mass ranges from approximately 0.5 GeV at m^^ = 10 GeV to 2.3 GeV at m^^ 
= 100 GeV. 
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4.2 General Search Results 

Figure |^ shows the di-photon mass versus the recoil mass for all candidate events passing the 
general search cuts for the year 2000 data only (where all the data were taken at £"0111 near 
206 GeV). The events at recoil masses near zero are expected from e''~e~ — > 77 background. 
This plot also shows no unexpected structure for the lower Ecm data. In the absence of an 
indication for signal, limits are placed on the production at £'cm~206 GeV of a massive state 
decaying into photon pairs. 

For the general search, the system recoiling from the di-photon system is not assumed to be 
a Z" and hence the branching fractions X ^ 77 and Y into topology A, B or C are not uniquely 
predicted. Here X is a scalar particle and Y is a scalar or vector particle. Furthermore, X and 
Y must be a short-lived particles so that they decay near the interaction point. In order to be 
independent of models we do not combine data from different Ecm and therefore we restrict this 
part of the analysis to the highest energy data, in the £"0111 range of 205 — 207 GeV; this represents 
200.0 pb-i of the 2000 data. We choose to present upper limits on a{e~^e — > XY) x B(X. — > 
77) X B{Y ^ ff) as a function of Mx. When presenting production upper limits as functions 
of Mx, we show the limit obtained for the value of My that gives the smallest efficiency in the 
region Mz — Mx < My < Ecm — Mx- The lower bounds on My are used because searches 
for di-photon resonances at LEPl ^ have already set good limits on the lower- mass 

phase space. Mx and My are also required to be above 10 GeV and below 200 GeV in order 
to allow the decay products to have sufficient energies and momenta to give reasonable search 
acceptances at Ecm = 206 GeV. For a scalar /vector hypothesis for X/Y, the efficiency is found 
to be the same to within 5% as that for a scalar/scalar hypothesis; the lower of these efficiencies 
is used in setting the limits. For the lepton search channel, the efficiency for Y — > r+r^ is used, 
as it turns out to have the lowest of the dilepton efficiencies. 

The event candidates from Ecm = 205 — 207 GeV in the general search are used to calculate 
95% CL upper limits on the number of events in 1 GeV [Mx,My] mass bins. The acceptances 
used at each 1 GeV mass bin are obtained by interpolation using a 4th-order polynomial fit to 
the acceptances simulated on a 5 GeV grid. For each [Mx,My] bin, the 95% CL upper limit on 
the number of signal events is computed using the frequentist method of reference [^8|, which 
takes into account the predicted Standard Model background. This statistical procedure also 
incorporates the di-photon mass resolution (typically less than 2 GeV for m^^<100 GeV); the 
limit procedure is valid for resonance states narrower than this resolution. The effect of the 10% 
systematic error for background modelling is incorporated in the statistical procedure, as is the 
4% uncertainty on signal. For these general search limits, an additional systematic uncertainty 
of 5% is added to the signal uncertainty to account for interpolation error in the efficiency grid 
(especially near kinematic limits) and for the differences in the acceptance calculations for the 
scalar versus vector nature of particle Y. 

Figure | shows the 95% CL upper hmits on cj(e+e- XY) x B{X 77) x B{Y ff). 
These results are valid independent of the nature of Y, provided it decays to two jets, a lepton 
pair, or missing energy, and has a width less than or equal to the experimental resolution. 
Limits of 25 — 60 fb are obtained over 10 < Mx < 180 GeV. The limits for the leptonic final 
state are stronger, except in the case Y couples exclusively to r+r^ (the final state with lowest 
acceptance) . 

4.3 Limits on h°Z° with 77 

The distribution of di-photon masses for the h^'Z" search candidates for the 1999 data {Ecm = 
192 - 202 GeV) and the 2000 data {Ecm = 200 - 209 GeV) is shown in Figure together with 



8 



the simulation of Standard Model backgrounds. The observation of 42 events is in reasonable 
agreement with the expected background of 51.9ib2.9 events. Because the h'^Z^ process has a 
production rate and branching fractions described by theory, the data taken at all LEP energies 
can be combined in this analysis. Figure ^3 shows the distribution of m-y^ for Ecm = 88 — 
209 GeV. This plot is restricted to ni^^ larger than 20 GeV because there is no background 
estimate for the low-m^^ LEPl data. The figure has no indication of a resonance, and the total 
of 124 candidates agrees with the predicted background of 135.2ibl0.8 events. 

Also indicated on Figure |3| is the hypothetical signal of a 100 GeV Higgs boson produced at 
Standard Model strength and decaying into di-photons with a branching fraction of 18% - the 
fraction predicted in the "benchmark fermiophobic model" calculated by simply turning off the 
Higgs-fermion coupling. In reality, the fermiophobic Higgs photon branching ratio depends on 
parameters and details of fermiophobic 2HDM models Q,^, so this benchmark is simply a guide 
to the broad interpretation of the data. Here we use the HDECAY package to calculate the 
modified photonic branching fractions. 

The events passing all h^Z^ cuts are used to set an upper limit on the di-photon branching 
ratio for a particle produced in association with a Z^ and having the Standard Model Higgs boson 
production rate. As described in the previous section, the frequentist method of reference Q| 
is used to determine the 95% confidence level upper limit on possible signal events at each di- 
photon mass. Figure Q shows the 95% CL upper limit for the di-photon branching ratio obtained 
by combining the candidate events in 1999 and 2000 data described in this paper with those from 
OPAL searches at ^ = 88 - 189 GeV |8|,|,|ll|, where the Standard Model h°Z° production 
cross-section is assumed at each centre-of-mass energy. Higgs bosons produced at Standard 
Model rate and decaying exclusively to di-photons are ruled out at the 95% confidence level over 
the mass range 20 — 117 GeV. Figure ^ also shows the h*^ — > 77 branching ratio computed using 
HDECAY with the fermionic couplings switched off; the photonic branching fraction falls as the 
W+W~ and Z'^Z'^ channels become kinematically favourable. A 95% CL lower mass limit for 
the benchmark fermiophobic Higgs bosons is set at 105.5 GeV, where the predicted branching 
ratio crosses the upper-limit curve. The median limit one would expect to obtain in an ensemble 
of experiments in the absence of a signal is 106.4 GeV. The benchmark fermiophobic branching 
ratios can also be calculated using the HZHA3 [29| generator. HZHA3 produces slightly higher 
di-photon branching fractions than does HDECAY. The lower mass limit on the benchmark 
fermiophobic Higgs boson calculated with HZHA3 is 106.3 GeV. 



5 Conclusions 

A search for the production of Higgs bosons and other new particles of width no larger than 
the experimental resolution and decaying to photon pairs has been performed using e"'"e~ 
annihilation data with Ecm = 192 — 209 GeV combined with 88 — 189 GeV data from previous 
OPAL searches. Model independent upper limits are obtained for £'cm~206 GeV on a{e'^e~ — > 
XY) X 5(X ^ 77) X 5(Y ^ ff), where hmits of 25-60 fb are obtained over 10 < Mx < 180 GeV, 
for 10 < My < 200 GeV and Mx + My > Mz- The limits are valid for Y either a scalar or 
vector particle, provided that the Y decays to a fermion pair (interpreted as two jets, a lepton 
pair, or missing energy). 

Using OPAL data from all LEP centre-of-mass energies, model-specific limits are placed on 
B{h^ — > 77) up to a Higgs boson mass of 117 GeV, provided the Higgs particle is produced 
via e'''e~ — > h'^Z'^ at the Standard Model rate. A lower mass bound of 105.5 GeV is set at the 
95% confidence level for benchmark fermiophobic Higgs bosons. Similar lower mass limits on 
benchmark fermiophobic Higgs bosons have been obtained by the other LEP experiments |12-Q. 
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Run (year) 


Integrated Luminosity (pb ^) 


^cm (GeV) 


1990-95 


173.00 


88 - 94 


1995 


5.41 


130 - 140 


1996 


10.32 


172.3 


1996 


10.04 


161.3 


1997 


57.73 


182.6 


1998 


182.61 


188.6 


1999 


28.90 


191.6 


1999 


74.79 


195.6 


1999 


77.21 


199.6 


1 non 
lyyy 


OD.Uo 


ZUi.O 


2000 


0.82 


200 - 202 


2000 


2.62 


202 - 204 


2000 


76.81 


204 - 206 


2000 


123.26 


206 - 208 


2000 


7.54 


208 - 210 



Table 1: Summary of all data sets used in the searches. Results using data from 1990 — 1998 
were reported in earlier publications [^Q. 
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Cut 


Data 


EBkgd 


(7/Z)* 


4f 




£100 (%) 


1999 Hadronic Channel 


(Al) 


10645 


10695.4 


7535.0 


3160.2 




70 


(A2) 


62 


56.1 


53.6 


2.6 




62 


(A3) 


48 


44.8 


42.5 


2.3 




61 


(A4) 


29 


22.0 ± 1.8 


19.8 


2.2 




61 


(A5) 


15 


10.9 ±0.9 


10.9 


0.0 




60 


2000 Hadronic Channel 


(Al) 


9371 


9152.4 


6096.0 


3056.2 




71 


(A2) 


52 


46.5 


43.5 


3.0 




60 


(A3) 


39 


38.3 


36.0 


2.3 




60 


(A4) 


18 


17.8 ±1.6 


15.7 


2.1 




58 


(A5) 


7 


7.9 ±0.7 


7.8 


0.1 




57 



Cut 


Data 


EBkgd 


e+e 


r+r- 




77 


e+e^ff 


£100 (%) 


1999 Leptonic Channel 


(Bl) 


41947 


39060.6 


37245.9 


715.6 


50.5 


314.6 


734.0 


82 


(B2) 


167 


188.6 


72.4 


11.3 


8.0 


95.5 


1.5 


69 


(B3) 


155 


178.8 


67.3 


10.4 


7.4 


92.4 


1.2 


63 


(B4) 


23 


31.4 ±5.8 


18.8 


5.0 


6.9 


0.4 


0.3 


50 


(B5) 


5 


9.2 ± 1.7 


4.6 


1.7 


2.9 


0.0 


0.0 


48 


2000 Leptonic Channel 


(Bl) 


37432 


33928.9 


32329.1 


607.8 


43.7 


273.9 


674.4 


82 


(B2) 


138 


146.6 


57.6 


10.2 


6.8 


70.9 


1.0 


71 


(B3) 


123 


141.7 


55.4 


9.2 


6.4 


69.7 


1.0 


67 


(B4) 


28 


23.7 ±4.8 


13.0 


4.4 


5.9 


0.2 


0.3 


61 


(B5) 


11 


9.8 ±2.0 


5.1 


1.9 


2.6 


0.0 


0.0 


52 




Cut 


Data 


EBkgd 




77 


e+e 


£+£- 


e+e-ff 


£100 (%) 


1999 Missing Energy Channel 


(CI) 


224989 


129713.4 


50.5 


3337.7 


124769.9 


157.4 


1397.9 


88 


(C2) 


377 


336.3 


13.0 


276.4 


45.1 


1.1 


0.7 


75 


(C3) 


71 


73.1 


12.2 


31.5 


27.8 


1.0 


0.5 


69 


(C4) 


33 


42.4 


12.1 


29.2 


0.9 


0.0 


0.2 


69 


(C5) 


8 


12.9 ±0.5 


11.4 


0.3 


0.9 


0.0 


0.2 


67 


(C6) 


3 


7.6 ±0.3 


7.5 


0.0 


0.0 


0.0 


0.1 


66 


2000 Missing Energy Channel 


(CI) 


202649 


113268.8 


43.9 


2737.6 


109070.9 


136.8 


1279.7 


87 


(C2) 


345 


315.1 


11.5 


267.5 


34.6 


0.9 


0.5 


75 


(C3) 


66 


62.1 


11.0 


30.1 


19.7 


0.8 


0.3 


70 


(C4) 


34 


39.0 


11.0 


27.8 


0.0 


0.0 


0.2 


69 


(C5) 


6 


10.5 ±0.5 


10.2 


0.1 


0.0 


0.0 


0.1 


68 


(C6) 


1 


6.5 ±0.3 


6.5 


0.0 


0.0 


0.0 


0.0 


65 



Table 2: The 77 searches: events from the 1999 and 2000 data remaining after the 

indicated cumulative cuts. Hadronic channel cuts (A), leptonic channel cuts (B), and missing 
energy channel cuts (C) are described in Sections 3.1, 3.2, and 3.3, respectively. The entries 
for A5, B5, C6 are for the Mrecoii cut for the h'^Z'^ search; they are not applied in the general 
search. The uncertainty shown indicates statistical error only; the systematic uncertainty is 
approximately 10%. In (A), the components from (7/Z)* and four-fermion ("4f' ) final states 
are shown. In (B), the components from Bhabha scattering (e^e"), r-pair, /i-pair, 77 and 
e+e-ff final states are shown. In (C), the components from h'V"y^, 77, e+e--pair, lepton pair 
{i = /u, r) production and e'^'e-ff final states are shown. The efficiencies for detection of a 
100 GeV Higgs boson are shown in the last column. 
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25 50 75 100 125 150 175 200 



(GeV) 

Figure 1: Distribution of mass recoiling against the di-photon system, Mrecoib versus di-photon 
invariant mass, m-y-y, for events passing the general search cuts on the Y2000 data. The different 
search channels are as indicated. The diagonal line denotes the kinematic limit. Dashed lines 
and arrows indicate the events accepted for the h'^Z'^ search. 
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Figure 2: For data taken at i?cm~206 GeV: 95% confidence level upper limit on a{e~^e~ — >■ 
XY) X i?(X ^ 77) X B{Y ^ ff) for the cases where Y decays hadronically (solid line), Y decays 
into a charged lepton pair (dashed line), and Y decays invisibly (dotted line). The limits for 
each Mx assume the smallest efficiency as a function of My such that 10 < My < 200 GeV and 
that Mx + My > Mz- 
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Figure 3: Distribution of mass of the two highest-energy photons in the h'^Z'^ search after 
apphcation of all selection criteria, a) shows the data taken in 1999 and 2000 only, while b) 
shows the data taken at all energies. Data are shown as points with error bars. Background 
simulation is shown as a histogram showing the contributions from the hadronic, charged lepton 
and missing energy channels as denoted. The expected signal for a 100 GeV fermiophobic Higgs 
boson produced with Standard Model strength is also shown. 



16 




Figure 4: 95% confidence level upper limit on the branching fraction B{h^ 77) for a Standard 
Model Higgs boson production rate. The shaded region, obtained with all LEP energies, is 
excluded. The dotted line is the predicted B{h^ —>■ 77) assuming B{h^ — > ff) = 0. The 
intersection of the dotted line with the exclusion curve gives a lower limit of 105.5 GeV for the 
fermiophobic Higgs model. The median expected limit and the ±20" range of expected limits 
are indicated by the dashed lines. 
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